A1 = readmatrix('NiI2_pol.csv')

x = (A1(:,1)*2+5)*2*pi/360
y = abs(A1(:,2))
z = A1(:,3)
y= y./z

x1 = (A1(1:20,1)*2+5)*2*pi/360
y1 = abs(A1(1:20,2))
z1 = A1(1:20,3)
y1= y1./z1

x2 = (A1(20:66,1)*2+5)*2*pi/360
y2 = abs(A1(20:66,2))
z2 = A1(20:66,3)
y2= y2./z2

x3 = (A1(66:110,1)*2+5)*2*pi/360
y3 = abs(A1(66:110,2))
z3 = A1(66:110,3)
y3= y3./z3

x4 = (A1(110:157,1)*2+5)*2*pi/360
y4 = abs(A1(110:157,2))
z4 = A1(110:157,3)
y4= y4./z4

x5 = (A1(157:181,1)*2+5)*2*pi/360
y5 = abs(A1(157:181,2))
z5 = A1(157:181,3)
y5= y5./z5

ax = polaraxes
polarplot(x1,y1,LineWidth=2,Color=[0 0 0.8])
hold on
polarplot(x2,y2,LineWidth=2,Color=[0.8 0 0])
polarplot(x3,y3,LineWidth=2,Color=[0 0 0.8])
polarplot(x4,y4,LineWidth=2,Color=[0.8 0 0])
polarplot(x5,y5,LineWidth=2,Color=[0 0 0.8])
polarscatter(x,y,5,'MarkerEdgeColor',[0 0 0])
ax.ThetaZeroLocation = 'top'
ax.ThetaDir = 'clockwise'
hold off

set(gcf,'Position',[0, 0, 200, 200])